Scattering of Carriers by Charged Dislocations in Semiconductors 
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The scattering of carriers by charged dislocations in semiconductors is examined within the frame- 
work of the Bohzmann transport theory. The ratios of quantum and transport scattering times are 
evaluated after averaging over the anisotropy in the relaxation time. A new approximate expression 
for the carrier mobility is proposed and the value of the Hall scattering factor is computed. Change 
in the resistivity when the dislocations are tilted with respect to the plane of transport is deter- 
mined. Finally an expression for the relaxation time is derived when the dislocations are located 
within the sample with a uniform angular distribution. 

PACS numbers: 72.10.Fk, 72.20.Dp 
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I. INTRODUCTION 

Epitaxial growth of thin semiconductor films on sub- 
strates which have a large lattice constant mismatch re- 
sults in the films being strained. Depending on the 
growth conditions and the films' thickness, this strain 
can either partially or fully relax through a formation of 
various possible kinds of lattice defects. Among these 
defects edge dislocations are prominent and have a pro- 
nounced effect on the mobility of carriers. While the the- 
ory for charged dislocation scattering was first formulated 
to explain the low temperature mobility of plastically de- 
formed semiconductors jH, interest in dislocation scatter- 
ing has revived in context of GaN In- 
deed it is important in all epitaxially grown materials 
on mismatched substrates, as well as bulk crystals whose 
growth techniques have not yet been mastered . 

An edge dislocation is a row of dangling bonds formed 
by an abruptly terminated plane somewhere inside the 
crystal 2J . This local departure from tetragonal coordina- 
tion produces acceptor states in the energy gap, forming 
one dimensional lines of charge. The effective screened 
electrostatic potential energy, U{x±) is thus cylindrically 
symmetric if the extent of the edge dislocation is taken 
to be infinite 113, El. 
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where Q is the charge per unit length, Kq is the modified 
zeroth order Bessel function of the second kind, e — eoCr 
is the dielectric constant, A is the screening length and 
is the distance from the dislocation line in a perpen- 
dicular plane, r = f{x_^^,9,z). These one dimensional 
lines of charge have detrimental effects on the transport 
properties of charge carriers. 



II. ISSUES ADDRESSED 

(i) The scattering potential, due to its cylindrical symme- 
try, is highly anisotropic and the validity of the extension 
of the relaxation time approach to this problem has been 
questioned[lJl[I3. 

(ii) Podor's [l| expression for the relaxation time 
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is apparently physically incorrect. wj_ is the component 
of electron velocity perpendicular to the dislocation axis 
and is the number of dislocations per unit area, all 
assumed to be parallel and independent. Only the per- 
pendicular component of the impinging electron's veloc- 
ity contributes to scattering and the component parallel 
to the dislocation is unaffected. Equation[21is finite when 
v±^ — > 0, whereas in this limit, r should diverge. 

(iii) The method of energy averaging employed by Podor 
has been questioned^. 

(iv) Due to (iii), the tensor nature of resistivity is not 
evident in the final expression. In particular if the dislo- 
cations are tilted at an angle with respect to the direc- 
tion perpendicular to the plane of transport, it is difficult 
to give anything better than a rough estimate^ in the 
present theory. Most often, the effect of dislocation ori- 
entation is disregarded and /i^ is replaced by a scalar 
number. 

(v) Quantum and classical scattering times were 
calculated E3 without averaging out the anisotropy in 
the problem. 

(vi) There are corrections to the measured Hall mobility 
due to the Hall scattering factor. 

(vii) In general, dislocations may not be all parallel. 

Points (i) and (ii) above were what led us to the re- 
examination of this problem. Nevertheless, a first prin- 
ciples treatment of the calculations shows that Podor's 
40-year old formula is indeed correct. The reasons for 
the apparent discrepancy are traced to the limits where 
the solution of the Boltzmann equation is invalid. This 
is discussed in the next section. 
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III. THEORETICAL FORMULATION 

We start from the Boltzmann equation within the lin- 
ear response regime (3 • Then upto first order in elec- 
tric field, the perturbed distribution function may sym- 
bolically be written as /k = fok + 0k 7 where 0k is de- 
viation from an equilibrium distribution in presence of 
a perturbing external electric field F. In absence of 
a thermal gradient and a magnetic field, the linearized 
Boltzmann equation for carriers described by spherical 
parabolic band reduces to 



m* ' dE 



E 

k' 



Wk\k[0k' - 0k] 
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Wk,k' is the transition rate between initial and final 
plane wave states, k and k', in presence of the scat- 
tering potential given by equation iQJ. For scattering 
from charged dislocations, the scattering rate is given by 
M^k.k' = <5(fc, - K)6{k - fc')5(|kl - k^l). g(|kl - k^l) 
is the part depending on only a function of in-plane mo- 
menta (shown below). Thus (a) collisions are elastic, (b) 
the components of the incident electron's momenta which 
are parallel and perpendicular to the dislocation line are 
separately conserved, (c) no electric field develops along 
the dislocation axis, i.e. F • k = Fj^ • kj^. This imme- 
diately implies that no relaxation time can be defined 
along the direction parallel to the dislocations' axis. In 
other words, for time independent electric field, there is 
no steady state solution to the Boltzmann equation if the 
collision term is zero. Nevertheless, one may physically 
argue that 1/tz = 0. The argument is clear within the 
variational formalism|l5j where one defines the sample 
resistivity in terms of the Joule-heat dissipated due to 
a finite current. With constraints (a)-(c) in mind, we 
shall choose a 0k which solves the Boltzmann's equation 
exactly. Ansatz: 



Here 9 is the angle between kj_, k'j_ which lie on a circle 
parallel to the xy-plane since is independently con- 
served. The wave vectors in the summation in equation 
(|7|) are three dimensional. 
Within the Born approximation 



WkM' = '^[Lz5k^.k'^''5{Ek-Ek') 
— i— / dxiC/(xi)e' 

JxLyLz J 



(ki-k'_^)-xj 
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Here Lz,Ly and are the crystal dimensions over 
which the plane wave electron states are normalized 
and the length of the 'infinite' dislocation has been lim- 
ited to the size of the crystal along the z direction. 
The energy conserving delta function, 5(Ek — Ek') = 
{dE/dk)-^S{k - k') = {nk/m*)-\k/k±)S{k± - k'^) due 
to Sk^. k'^ in the summation. Thus, as previously claimed, 
both the perpendicular and the parallel components of 
the electron momenta are separately conserved. Since, 
Ek' LxLyl {2'kY ^ dk'^, an overall factor of area re- 
mains in the denominator after the primed momenta have 
been integrated over. This simply means that the scat- 
tering due to a sing le charged dislocation is ineffective in 
a large sample, [l^ When there are many charged dislo- 
cations within this area which are all parallel, one can 
simply replace {L^Ly)^^ by Nd the dislocation density 
per unit area if the interference terms can be neglected. 
The Fourier transform in equation (jSJ can be approxi- 
mately written as Q 



C/(|kl 



QeA^ 



e(l + |k'^-ki|2A2) 



(9) 



IV. RESULTS 



A. Transport Lifetime 



m" dE 
Substituting 0k in equation© yields 



(4) 



Ok ^ 



dE 



dk 



Ok' 



dE 



r(fc',fcl)k'j 



Ok 



dE 



T{k, fc±)kj 
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From energy and perpendicular momentum conservation, 

(6) 



^^°'='-r(fc',/fcl) = ^r(fc,fc^) 



dE 



dE 



Thus the linearized Boltzmann equation is exactly solved 
if 



(7) 



From equationQ, the relaxation time in the direction 
perpendicular to the dislocation axis is 



T^{k, k±_) 



Q'^e^X'^m*Nd 



(10) 

This is exactly what Podor had derived (equation Q) 
and k± — implies a finite tj_ even after our rederiva- 
tion. While in three dimensions an electron with fc = is 
unphysical (there is no associated phase space), an elec- 
tron with k± = and kz ^ corresponds to a physical 
situation. The inconsistency in the final formula results 
from the breakdown of the validity of the assumed solu- 
tion, 0k = for k± = in equation Q). This condition 
is outside the scope of the present scheme of the solution, 
which is otherwise consistent. 

The anisotropy in t necessitates a further angular av- 
eraging for a comparison with any physical quantity as- 
sociated with a measurement which involves a thermody- 
namic distribution of electrons. This transport scattering 
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time is directly connected to mobility, fj, = (e/m)((r)), 
where (()) denote an energy average, (see below) over 
a distribution function of appropriate degeneracy. In a 
fully degenerate system, using equation , this simpli- 
fies to {{Ttr)) = (3/4) /J" sin^ 6Itj_ d9. 



B. Quantum Scattering Time 



A quantum scattering time, T'^{k±) is, by definition, 
equation ||7J), but without the (1 — cos 6) factor and may 
be calculated similarly. This was recently done by Jena 
and Mishrafn. 



[1 + (2fciA)2]3/2 



(11) 



However, the angular dependence of must also to 
be averaged out. The meaningful quantity is (I/t'^) = 

{2/tt) J^^'^[T'^{9)]~^d6 and is often connected to the fi- 
nite amplitude and width of the Shubnikov-de Haas or 
de Haas-van Alphen oscillations. The quantum scatter- 
ing time may be looked upon as an effective 'Dingle' tem- 
perature, Td ~ (?i/27rfcB)(l/r9). 

Nevertheless, while comparing Shubnikov amplitudes, 
the scattering rates are better calculated between Lan- 
dau wave functions and with a density of states at the 
Fermi level modified by the magnetic field, as was done 
long back by Vinokur[l7l| for the essentially the same 
problem. Furthermore, literature on the connection be- 
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FIG. 1; The ratio of dislocation scattering limited transport 
and quantum scattering times for a degenerate electron gas. 



tween scattering times for dislocations' strain field and 
de Haas- van Alphen oscillation amplitudes in metals was 



a subject of lively debate sometime back. Many par- 
allel interpretations for level broadening [l^l have been 
suggested. Some semiclassical arguments even favour a 
small angle cutoff. This fact may be particularly impor- 
tant in two dimensions where it could rescue the quan- 
tum scattering time from a divergence in a simple 
and physically meaningful way, the small angle cutoff 9c 
(in radians) being invers ely pro portional to the Landau 
level index n, 9c ^ n/2n 

Despite the preceding remarks, the concept of a quan- 
tum scattering time finds a widespread use in literature 
(for example, references poll2lLl2^'). Therefore we have 
plotted the suitably defined ratio / t"^) {{nr)) of the 
transport and quantum scattering times for a three di- 
mensional degenerate carrier gas in figure ^ The graph 
is plotted as a function of the dimensionless parame- 
ter, kp/qxF- Qtf is the simple wave vector independent 
Thomas-Fermi screening function. The largeness of this 
ratio is often regarded as a measure of 'anisotropy' of 
scattering. 5^. The real space anisotropy of the dis- 
location potential is different from the anisotropy in its 
Fourier transform, which is more a measure of the ef- 
fective range of the potential. An additional averaging 
causes the transport to quantum scattering times ratio 
to be larger than what was calculated by in reference 

El. 



C. Mobility 

In calculating mobility, the averaging procedure em- 
ployed by Podor has been called 'unspecified' f3| and 
hence it is worked out below. For dislocations along the 
z-axis, the current and electric field directions coincide as 
long as the measurement is done in the xy-plane. Then, 
jx ^ ne((vx)) and {(vx)) = M^Fx where 



{{Vx)) 



Or 



Or 



1*2 



Ek(/o + </'k) Ek/Ofe 



-^]r^(k,k^)d'k 



(27r)3 



(12) 



(13) 



d9 sin^ ( 



dkk 



4/ dfo^r 

dE I L 



(2A:Asin6') 



21 3/2 
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The integrals must now be evaluated numerically. Equa- 
tion (|13|l has the unpleasant feature of depending very 
strongly on screening length and thus at low tempera- 
tures turns out to be dependent on the model used for 
the temperature dependent of carrier concentration and 
screening. A simple analytic expression can be estimated 
by interpolating the two integrals (/ d9 and / dk) be- 
tween the two extremes cases, when the first term is much 
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smaller and when it is much larger than the second term 
in square brackets in equation (|13|1 . This is significantly 
better than Podor's high temperature approximation [ij 
(fcAsin6' ;s> 1). The relative percentage errors are plotted 
in figure [3 as a function of the dimensionless parameter 
8m i^bTX ^ seen that this approximation of the 

integral never deviates from the numerically calculated 
exact answer by more than 5%. 



jx = f^xxEx + cTxyEy, the Hall scattering factor rn is 
defined as 



rn = ne 



^ xy 



(17) 



where the off-diagonal conductivity, Uxy^ for carriers with 
parabolic energy dispersion which are distributed along 
isotropic constant energy surfaces is 
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in our formula and Podor's approximation wi§i respect to the 
exact expression evaluated numerically. The graphs are plot- 



ted as a function of dimensionless parameter 

Assuming that the electrons are distributed according 
to Maxwell-Boltzmann distribution, 



47ri/2?i3^2 



rl/3 



3/2 
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and when the carrier gas is fully degenerate 

(4/3)'/^ (57r/16)'/' Ak^pX^TP 

(16) 
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D. Hall Factor 

In most experiments it is not the drift but the Hall 
mobility which is measured. Under the assumption that 
the scattering rate does not alter in presence of a mag- 
netic field, B and when the magnetic field is aligned with 
the dislocations' axis, only the in-plane relaxation time 
comes into the picture. Using the same line of argu- 
ments, it is easy to again establish its existence for ar- 
bitrarily strong non-quantizing magnetic fields. Then, if 



e^B f (fik 2^/0 f dE'^'^ 



'xy — 2 
h 7 



47r3 ^dE \dk 



1 + 

\ m* J 



(18) 



From equations (|13|l . (|18|l and H17|) the Hall scattering 
factor for nondegenerate carriers at high temperatures 
(i. e. 8'"*fef > 1) approaches a value of 2.17, ob- 
tained by dropping the second term in square brackets 
in equation IjlSII . At lower temperatures, its value is de- 
pendent on the model of carrier density and screening 
but always smaller. The anisotropy in scattering makes 
the value higher than the Hall factor for ionized impurity 
scattering which is 1.93. 



E. Effect of Dislocation Tilt 

Assume that dislocations are all parallel, but now 
at a longitude (j) and latitude 9 with respect to the z- 
axis while the measurement is being done in the xy- 
plane. A unit vector along this dislocation axis is d = 
X sin 9 sin (j) + y sin 9 cos (j) + z cos 9. Because the electric 
field is developed only along the direction perpendicular 
to the dislocations' axis, Fj_ — pi± = p[j — (j ■d)d] which 
yields (with cos and sin abbreviated to c and s) 



P =P 



l~s^9s^(f> -s^9s(f)c4 
~s^9s(t)C(j) 1 - s^9c'^(, 
—c9s9s(l) —c9s9s(h 



-~c9s9sq 
-~c9s9s(f 
l-c^9 



(19) 
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Negative sign in the off diagonals indicates the direction 
of the electric field developed. 



F. Angular Distribution of Dislocations 

Let us next consider the case when all the dislocations 
are not parallel to each other but are instead distributed 
with a uniform distribution of angles. One can, of course, 
average over the angles(2Jj appearing in equation H19|l . 
This averaging over the angles in the rotated resistivity 
tensor amounts the use of Matthiessen's rule and will not 
change the temperature dependence of mobility. 

For a better approximation, we again start from the 
linearized Boltzmann equation, equation (jS)). In the 
present case, the relaxation time must be isotropic and 
therefore let the ansatz for the distribution function be 



m 



he dfi 



Ok 



m dEh 



Tisoik)k ■ F 



(20) 
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We shall further assume incoherent scattering such the 
scattering rates due to different dislocation lines add. If 
the scattering rate due to an dislocation is j^., , then 

the total rate is J2i ^k,k'- 

Without loss of generality, one can choose the electron 
wave vector k to be along the z-axis, k — kz. If the axis 
of the i*'* dislocation, d* is at an angle {6, (f)) with respect 
to the z-axis, then the unit vector along the dislocation 
axis is given by = smO sm(j)x + sin0cos(/)y + cos^z. 
The component of the wave vector perpendicular to the 
dislocation axis is given by 



Or 



ki = k - (k-d')d' 



k [(1 — cos^ 9) z — cos 9 sin 9 cos (j) y 



— cos 6* sin sin x] (21 
Substituting back in the Boltzmann equation, we get 

1 



k* 



T(fc,/Cl) 



(22) 



i 

Converting the sum into an integral, 

^.fc. = Wfc)^F./df^kl-^ (23) 
Here Na is the dislocation density per unit solid angle. 



Since the averaging over the dislocation orientations is 
equivalent to an averaging over the electron wave vectors, 
the expression for the relaxation time becomes 



1 



Tiso {k) 



47r 



d9- 



(1 



cos 



[l + fc2A2 sin^ 6l]3/2 
(24) 

where the integral has been performed and we have 
noted that J^^ sin (j}d(j) = J^^ cos 4>d4> = 0. From here on, 
it is straightforward to calculate the isotropic mobility, 
although it must be done numerically psf . 



V. SUMMARY 

Within the framework of the conventional transport 
theory, we have shown that a relaxation time can be 
defined for scattering of carriers by charged disloca- 
tions. Difference between quantum and classical scat- 
tering times was discussed and it was pointed out that 
the anisotropy necessitates an appropriate angular av- 
eraging. A new approximate formula for mobility was 
derived and it was shown be within 5% of the exact re- 
sult at all temperatures. The value of the Hall scattering 
factor and the effect of dislocation tilt on resistivity was 
determined. Finally we derived a new expression for the 
relaxation time when the angular orientation of disloca- 
tions is isotropic. 
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